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ABSTRACT 

We present the discovery of five millisecond pulsars (MSPs) from the PALFA Galactic plane survey 
using Arecibo. Four of these (PSRs J0557+1551, ji850+0244, JI902+0300, and JI943+2210) are 
binary pulsars whose companions are likely white dwarfs, and one (PSR J1905+0453) is isolated. 

Phase-coherent timing solutions, ranging from ~ 1 to ~ 3 years in length, and based on observations 
from the Jodrell Bank and Arecibo telescopes, provide precise determinations of spin, orbital, and 
astrometric parameters. All five pulsars have large dispersion measures (> 100pccm“^, within the 
top 20% of all known Galactic field MSPs) and are faint (1.4 GHz flux density ;$0.ImJy, within 
the faintest 5% of all known Galactic field MSPs), illustrating PALFA’s ability to find increasingly 
faint, distant MSPs in the Galactic plane. In particular, PSR J1850-I-0244 has a dispersion measure 
of SdOpccm”®^, the highest of all known MSPs. Such distant, faint MSPs are important input for 
accurately modeling the total Galactic MSP population. 

Keywords: pulsars: general — pulsars: individual (PSR J0557-I-1551, PSR J1850-I-0244, 

PSR JI902-h0300, PSR J1905-f0453, PSR J1943-f22I0) 


1. INTRODUCTION 

Millisecond pulsars (MSPs) are neutron stars that 
have been spun up by accreti on to ext r emely 
fast spin periods (P < 20 ms ) (jAloar et al.l 11981^ 
Radhakrishnan fc SrinivasanllI98^ . The first discovered 
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MSP , PSR B1937+21, was found in 1982 ([Backer et alJ 
Il982f) and since then ^160 Galactic MSPs have been dis¬ 
covered (~130 MSPs have been found in globular clus¬ 
ters) in different pulsar surveys (f or a review of pas t 
and ongoing searches for MSPs see iStovall et al.ir2013[) . 
MSPs are of great scientific value, giving us insight into 
many different branches of fundamental physics and as¬ 
trophysics, as we elaborate below. 

Most MSPs are found in binary systems with a low- 
mass, degenerate companion. These companions are 
responsible for spinning-up such pulsars to millisec¬ 
ond rotational periods by transferring mass and an¬ 
gular momentum onto them. The discovery of new 
MSP binaries can test (and in some cases confirm) 
mode l s of binary evolut i on (e . g. iPhinnev fc Kulkarnil 
1994 iTauris fc Savoniid 1199^ lArchibald et al.l 120091: 
Paoitto et al.l I2013D md sometimes challenge them, 
as in the case of PSR J1903-h0327 ([Champion et alJ 
12008(1 . which suggested that formation in a triple 
stella r system occurs in some cases ([Freire et al.l 
1201 ll: iPortegies Zwart et al.l 12011(1 . This was also di¬ 
rectly demonstrated in t he recently discove red PSR 
J0337-I-1715 triple system ([Ransom et al.l[^014[) . On the 
other hand, isolated MSPs exis t as well, and their origin 
is an important open question (iBai lyii fc G rindlavlll99rt 
iLevinson fc EichlCTll991l : iBailes et al.ll2nil[L 

The continued discovery of MSPs is also motivated 
by their potential for aiding the detection of gravita¬ 
tional waves. High-precision, lon g-term timing of MSPs 
forming a ‘pulsar timing array’ ((Foster fc Backedll990H 
could allow the detection of gravitational waves in the 
nanohertz frequency range particularity in the form of a 
stochastic background that is due to super-massive black- 
hole mergers at high redshifts (fvnn Haasteren et al.ll2011l: 

Yardlev et al.l[2MlLlDemorest et al.ll2013HShannon et all 

2013f) . The discovery of more MSPs that can be 
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timed with very high precision is crucial for this ef¬ 
fort to succeed. Some MSP binaries also enable tests 
of General Rela tivity and alternate theories of relativis¬ 
tic gravity ('e.g iGonzalez et al.l[20Tlt iFreire et al.ll2012t 
lAntoniadis et al.ll201,‘lli as well as measure precise compo¬ 
nent m asses that can test models of neutron-star interi¬ 
ors (e.g.lDemorest et al.ll20l^ iLattimer fc Prakashll201Cll : 
iBednarek et al.ll2012H . 

PALFA is a pulsar survey using the 7-beam Arecibo 
L-Band Feed Array (ALFA) on the Arecibo Observatory 
William E. Gordon 305-m Telescope. The survey’s centre 
frequency is 1.4 GHz, and the planned survey coverage 
is the region of the Galactic plane (|6| < 5°) accessible 
from the Arecibo telescope (32° < I < 77° and 168° 

I < 214°). The survey began in 2004 and is ongoing 
(jGordes et al.ll200^ |L azaru ij boLsn . 

PALFA provides a unique window for the discovery of 
MSPs. Its high sensitivity allows observations with short 
integration times which mitigate the effects of Doppler 
smear ing of the pulsed signal in th e case of a tight binary 
orbit ((Johnston fc Kulkarnil[T9M^ . The increased sensi¬ 
tivity and relatively high observing frequency - which re¬ 
duces the deleterious effects of dispersion-smearing and 
multi-path scattering - compared to other surveys also 
allow PALFA to find MSPs at much further distances 
and higher dispersion measures (DMs). Figure [T] shows 
the unique parameter space probed by PALFA. 

The discovery of distant, highly dispersed MSPs will 
allow a more complete characterization of the Galactic 
MSP population. Gurrent models of the Galactic MSP 
population are based p rimarily on the Park es Multibeam 
Pulsar Survey ((PMPS: iLorimer et al.lf2f)?)^ which is sig¬ 
nificantly less sensitive than PALFA and therefore bi¬ 
ased towards nearby sources. The increased sensitivity of 
PALFA not only allows higher DMs, and thus distances, 
to be probed, but allows discoveries of nearby pulsars to 
a lower luminosity. The MSPs discovered in PALFA will 
therefore provide valuable input to models of the Galac¬ 
tic MSP population - especially the population that is 
confined to the Galactic plane. 

High-DM MSPs also provide unique probes into the 
properties of the interstellar medium (ISM). As pul¬ 
sar signals propagate through the ISM, their signals are 
both dispersed (frequency-dependant delay), and scat¬ 
tered (multi-path propogation). These effects are time 
variable because of the pulsar’s proper motion and thus 
our changing line-of-sight through the ISM. Detailed 
dispersion and scattering measurements of signals from 
high-DM MSPs can be used to compare the observed 
effects to the predicted, thus improving our models of 
scattering and time-dependant DM variations. This is 
important for the detection of gravitional waves, as cor¬ 
recting the effects of DM variations and scattering can 
greatly increase the prec ision with which we time some 
MSPs (iKeith et al.ll20T^ . 

PALFA has discovered 18 Galactic plane MSPs 
to date, including seven t hat have been pub¬ 
lished in previous works (iGhampion et al.l I2008I 
iDeneva et alJ l2012t iCrawford et al.l I2012D . These in¬ 
clude: PSR J1950-I-2414, an eccentric MSP binary 
that may have formed vi a accretion induced c ollapse 
(Knipsel et al. in prep; IFreire fc TanrisI 1201411 : and 
two MSPs for which Shapiro delay measurements were 
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Figure 1. Period vs dispersion measure fo r all known Galactic 
MSPs (cf. Figure 5 of ICrawford et aLl 120121) . PALFA provides a 
unique window at high DMs and f a st spi n periods. The PALFA 
pulsars are fro nilG^ niDion et aLI l|2Q08l) : IDeneva et al.l H2012 L 
ICrawford et al.l l|2012IL The ‘ATNF’ MSPs are for all Galactic 
MSPs listed in the ATNF pulsar catalog l|Manchester et aLll2005l‘) ‘^. 


made, PSR J190 3-I-0327 (IGhampion et al.l l2008fl and 
PSR J1949-f3106 (IDeneva et al.ll2012D . 

In this paper we present five PALFA MSP discov¬ 
eries: PSRs J0557-bl551, J1850-b0244, J1902-b0300, 
J1905-I-0453, and J1943-I-2210. The remaining six un¬ 
published PALFA-discovered MSPs will be presented in 
forthcoming papers along with future discoveries. Sec¬ 
tion [5] describes the discoveries as well as the follow-up 
timing observations and analysis. In Section |3] we dis¬ 
cuss the results of our timing analysis and explore the 
interesting properties of these newly discovered pulsars. 
In Section 4, we summarize our conclusions. 

2. OBSERVATIONS & ANALYSIS 
2.1. Discovery 

The five MSPs presented here were discovered in 
PALEA survey observations taken between 2011 Septem¬ 
ber and 2013 June. The PALFA survey currently uses 
the Mock spectrometers^® as a backend for the 7-beam 
ALFA receiver at a central frequency of 1375 MHz. The 
Mock backend provides 322.617 MHz of bandwidth and 
a time resolution of 65.476 ^s for each beam of ALFA. 
PALFA has survey regions towards both the inner and 
outer Galactic plane. The outer-galaxy pointings are 
provided by our commensal partners, the ALFA Zone- 
of-Avoidance Survey, who are mapping extragalactic 
neutral hydrogen ern ission behind the Galactic plane 
((Henning et al.l I2O1O0 . A typical inner-galaxy survey 
pointing is 268 s and the outer-galaxy pointings are typ¬ 
ically 180 s in length. Four of the five pulsars presented 
here were discovered using the above setup with the Mock 
backend in the inner galaxy. PSR J0557-I-1551 was dis¬ 
covered with the Mock backend in an outer-galaxy point¬ 
ing. The discovery observation dates for each pulsar are 
listed in Table [T] 

PALFA survey observations are processed by several 
pipelines. O ne of these is based on the PRESTO soft¬ 
ware package ((Ransom et al.ll20(i^ ^^ and is described in 

http://www.atnf.csiro.au/research/pulsar/psrcat/ 
http://www.naic.edu/-astro/mock.shtml 
http://www.cv.nrao.edu/~ sransom/presto/ 
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Lazarus et al. (in prep). In brief, the PRESTQ-based 
pipeline operates on raw data automatically downloaded 
from a PALFA-administered archive at the Cornell Cen¬ 
ter for Advanced Computing and is processed on Guil- 
limin, a Compute Canada/Calcul Quebec facility oper¬ 
ated by McGill University. Some processing was also 
performed on dedicated computer clusters at the Univer¬ 
sity of British Golumbia and McGill University. Searches 
are performed both in the Fourier domain, for periodic 
sources, and in the time domain, for single pulses, us¬ 
ing standard PRESTO tools. PALFA observations are also 
searched by the Einstein@Home pulsar search pipeline 
which uses distributed volunt eer computing and is de¬ 
scribed in lAllen et all (120131) . Lastly, a ‘quicklook’ 
pipeline (also PRESTO-based), which searches reduced- 
resolution data in near realtime at the telescope, is also 
employed. Four of the pulsars in this work were discov¬ 
ered using the full-resolution PRESTO-based pipeline and 
one, PSR J1943-I-2210, was discovered with the quicklook 
pipeline. 

2.2. Timing Observations 

Once confirmed as pulsars, we commenced timing ob¬ 
servations of the five new MSPs in order to determine 
their rotational ephemerides. Timing observations be¬ 
gan at the Lovell telescope at Jodrell Bank Observatory 
shortly after discovery and Arecibo observations began 
up to a year later. 

Follow-up observations at Jodrell Bank used a dual¬ 
polarization cryogenic receiver on the 76-m Lovell tele¬ 
scope, having a system equivalent flux density of 25 Jy. 
Data were processed by a digital filterbank which covered 
the frequency band between 1350 MHz and 1700 MHz 
with channels of 0.5-MHz bandwidth. Observations were 
typically made with a total duration of 10 - 60 minutes, 
depending upon the discovery signal-to-noise ratio. Data 
were folded at the nominal topocentric period of the pul¬ 
sar for subintegration times of 10 seconds. After “clean¬ 
ing” of any radio-frequency interference (RFI) using a 
median zapping algorithm, followed by visual inspection 
and manual removal of any remaining RFI, the profiles 
were incoherently dedispersed at the nominal value of the 
pulsar DM. Initial pulsar parameters were established by 
conducting local searches in period and DM about the 
nominal discovery values and finally summed over fre¬ 
quency and time to produce integrated profiles. 

The Arecibo observations were performed with the 
L-wide receiver which has a frequency range of 1150- 
1730 MHz and a system equivalent flux density of 3Jy. 
The Puerto Rico Ultimate Pulsar Processsing Instru¬ 
ment (PUPPI) backend was used to record the data. 
Initially, each pulsar was observed in incoherent search 
mode with PUPPI, where a full 800-MHz-wide spec¬ 
trum in 2048 channels is read out every 40.96 fxs. This 
produces a high data volume, so each observation was 
post-processed by downsampling to 128 channels while 
correcting for delays due to dispersive effects from the 
interstellar medium within each channel before sum¬ 
ming. Later observations were performed with PUPPI 
in coherent-dedispersion mode (hereafter referred to as 
coherent mode) where pulse prohles were folded in real¬ 
time and recorded in 10-s integrations with 512 frequency 
channels. Data from outside of the 1150-1750 MHz L- 
wide band were removed. Radio frequency interference 


was excised with both au tomatic algorithm s and manu¬ 
ally using the PSRCHIVE (jHotan et al.ll200^ ^^ tools paz 
and psrzap. Typically, each pulsar was observed for 10 
min per observing session, except for PSR J0557-I-1551 
which was observed for ~30 min-1 hour per session. 

2.3. Timing Analysis 

Initial timing solutions were obtained using solely Jo¬ 
drell Bank observations. For the binary MSPs, orbital 
solutions were found by measuring the frequency of the 
pulsar in each observation and fitting a model of orbital 
Doppler shifts for a circular orbit. These orbital solu¬ 
tions were used as a starting point in the iterative timing 
procedure. Pulse times-of-arrival (TOAs) were obtained 
after matching with a standard template constructed 
from the Jodell Bank observations and processed us¬ 
ing standard pulsar timing techniques with PSRTIME and 

TEMPO^^^ 

The Arecibo observations were processed using stan¬ 
dard PSRCHIVE tools. Using pat, TOAs for each of 
the early Arecibo observations were extracted by c ross- 
correlating in the Fourier domain (jTavloil 119921 ) the 
folded profile with a noise-free Gaussian template de¬ 
rived from the first Arecibo observation. TOA errors 
were derived using a Markov chain Monte Garlo method 
(FDM algorihm in pat). These TOAs were then fitted 
with timing solutions using TEMPO. Our timin g analysis 
used the JPL DE405 solar-system ephemeris ( Standisbl 
I1998D and the UTC(NIST) time standard. 

In order to obtain the final timing solution presented 
here a high signal-to-noise template for each pulsar was 
constructed. To construct the templates, the Arecibo 
subbanded search mode observations were folded with 
the best ephemeris. A high signal-to-noise profile was 
then constructed by summing the profiles together, 
weighting by the signal-to-noise ratio of each observation. 
The summed profiles were subsequently smoothed using 
the PSRCHIVE tool psrsmooth to create a final template 
profile. Eor PSR J0557-I-1551, the procedure outlined 
above was used to construct a template profile, but us¬ 
ing the coherent-mode data instead. Figure [2] shows the 
summed profiles and smoothed templates. The smoothed 
templates were used to extract TOAs that were used in 
the final TEMPO fit. 

For nearly circular binary systems, there is a high cor¬ 
relation between the longitude of periastron, w, and the 
epoch of periastron passage, Tq, in timing models. Since 
the four binary pulsars presented here are nearly cir cu- 
lar, we use the ELLl orbital model (|Lange et al.ll200il) in 
TEMPO. It is parametarized with ei = esinw, 62 = ecosw, 
and Tasc = R) — wPt,/27r, where e is the eccentricity, Tasc 
is the epoch of the ascending node, and Pb is the orbital 
period. This parameterization breaks the covariance but 
it is an approximation to first order in e and so is valid 
only when xe^ is much smaller than the TOA precision. 
This is true for all the binary MSPs we are timing here. 

The TOA uncertainties for each telescope and data- 
taking mode (Jodrell Bank, Arecibo incoherent-search 
mode, and Arecibo coherent mode) were increased by a 
scaling factor (EFAC) that yields a reduced equal to 
unity. This results in more conservative estimates for the 

http://psrchive.sourceforge.net/ 

http://tempo.sourceforge.n6t/ 



















4 


X 


X 

ptH 


X 

PlH 




Figure 2. 1.4-GHz pulse profiles for the MSPs presented here. One full period is shown in 256 bins for each pulsar. The profiles were 
constructed by phase aligning and summing the Arecibo observations. The heavy black line shows the summed profiles and the light red 
line shows the profiles after smoothing was applied. The smoothed profiles were used as template for the generation of TOAs. 


timing parameter uncertainties. This is standard prac¬ 
tice in pulsar timing and it is well known that formal 
TOA uncertainties are often underestimated. These scal¬ 
ing factors (EFACs) are listed in Table [1] along with our 
best-ht timing solutions and derived parameters. The 
residuals are shown in Figure [H The quoted uncertain¬ 
ties in [T] are 68% confidence limits as reported by TEMPO. 

To make a refined measurement of the DM for each pul¬ 
sar, we split the 600-MHz bandwidth of each PUPPI ob¬ 
servation into six frequency bands and derived six TOAs 
from each observation for each 100-MHz band. We took 
our best-fit solution for each pulsar and held all param¬ 


eters except for DM fixed. We then fit the TOAs using 
all frequency bands to measure a DM. DMs measured in 
this fashion are listed in Table [TJ 
We also ht for proper motions for PSRs J0557-I-1551, 
J1902-h0300, and J1905-h0453. We found that for all 
three sources the measured values were consistent with 
no proper motion in either right ascension or declination 
within 95.4% (2cr) conhdence intervals. The upper limits 
on the total proper motion are shown in Table [TJ The 
timing baselines for PSRs J1850-I-0244 and J1943-I-2210 
were not long enough to break degeneracies with other 
parameters and place meaningful limits on the proper 
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Figure 3. Timing residuals for the five MSPs. Gray points indicate Jodrell Bank TOAs, red points denote PUPPI search-mode 
TOAs and blue points represent PUPPI coherent-mode TOAs. The inset shows a zoom-in of the Arecibo coherent-mode observations for 
PSR J0557-I-1551. The x-axes of the inset are aligned with those of the main plot. Error bars are scaled so that the best-fit reduced is 
equal to unity. The scaling factors along with rms values for the residuals are presented in Table [J 


motion. 

2.4. Flux and Polarization Measurements 

The Arecibo coherent-mode observations were flux- 
and polarization-calibrated using observations of a 
noise diode and of bright quasars (J1413-I-1509 and 
B1442-I-10) obse rved for the NANOGrav collaboration 
(|Demorest et al.l[2013f) . The calibrations were performed 
using the PSRCHIVE tool pac and the SingleAxis model, 
which assumes that the two polarization receivers are 
perfectly orthogonal. These calibrated observations pro¬ 
vide a full set of Stokes parameters for each 10-s inte¬ 
grated pulse profile and so measurements of the flux and 
polarization characteristics are possible. Observations in 
this mode were made for four of the five pulsars, with 


PSR J1850-I-0244 being the exception due to its more re¬ 
cent discovery. The set of calibrated coherent-mode ob¬ 
servations were summed and Figure |4] shows the summed 
polarization profiles for each of the four pulsars. 

In order to measure and correct for the effect of Fara¬ 
day rotation on the polarization of the pulsars, rota¬ 
tion measures (RMs) were determined by summing the 
profiles in eight frequency sub-bands for all the cali¬ 
brated coherent-mode observations. A set of trial RMs 
was tested by appropriately rotating Stokes U and Q 
for each frequency channel and then measuring the to¬ 
tal linearly polarized flux, L. The RM and its uncer¬ 
tainty were then measured from the RM-L curve. Only 
PSR J0557-I-1551 provided a significant measurement 
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Figure 4. Polarization profiles for the four MSPs that had data taken with full Stokes parameters. In the bottom panel of each plot, the 
black, red, and blue lines show the total, linearly polarized, and circularly polarized flux, respectively. The top panels show the position 
angle of the linear polarization for profile bins with a linear polarization having a signal-to-noise greater than three. An RM could only be 
determined for PSR J0557+1551 and so its polarization has been corrected for Faraday rotation whereas the others have not. 



Phase 


Figure 5. Pulse profile of PSR J1850+0244 (black line) fitted 
with a wrapped exponentially modified gaussian (red line). The 
unscattered gaussian, i.e. the underlying profile of the pulsar is 
represented by the dashed grey line. See text for parameters of the 
fit. 

with RM = 26 ± lradm“^ where the uncertainty is 
statistical. This value may co ntain a contribut ion from 
ionospheric Faraday rotation. iHan et ahl (|2006f) find for 
RMs for 223 pulsars at 1.4 GHz that the ionospheric RM 
is typically between —1 and — 5radm“^. So we add in 
quadrature an uncertainty of ±3 rad m“^ to the measure¬ 
ment error to get RM = 26±3radm“^. The polarization 
profile of PSR J0557-I-1551 in Figure|4]has been corrected 


for Faraday rotation using the measured RM value. 

The calibrated observations were also used to measure 
the flux densities of the four pulsars. An off-pulse re¬ 
gion was manually selected and its mean was subtracted 
from the calibrated profiles. For PSRs J1902+0300 and 
J1943-I-2210, the phase range excluding ±15% of pulse 
phase from the peak was used as the off-pulse region. 
For J0557±1551, due to its narrower pulse, ±10% from 
the peak was excluded. For J1905±0453, 30% of phase 
preceeding the peak and 40% following the peak was ex¬ 
cluded. The mean flux density of the baseline-subtracted 
profile was then computed. Those values are presented 
in Table [TJ Note that the quoted uncertainties are only 
statistical and do not take into account any uncertainties 
from the calibration procedure. 

3. DISCUSSION 

We have reported the discovery of five MSPs from the 
PALFA survey. We presented phase-coherent timing so¬ 
lutions of the five pulsars using observations from the 
Arecibo and Lovell telescopes. One of the pulsars was 
found to be isolated, and the rest are in nearly circular 
binaries. The DMs of the MSPs are high (see Figured]), 
especially PSR JI850±0244 (DM = 540pccm“^), which 
is the highest DM MSP yet discovered. PALFA has to 
date discovered 18 MSPs. All but one, and all five pre¬ 
sented in this work, have DMs greater than 100pc cm“^. 
This range is in the top 15% of non-PALFA Galactic 
















7 



Orbital Period (days) 


10 


10 


^ 2 
^ 10 


C 1 ri^ - O n 8 ° □ ° 

0 10 ° nn ^ ° 

Ph o ° “ 8°^ ° "J1943+2210 

^ J1902+0390„ s° ' J0557 + 1551 

^ A o o o°° a °o 

13 lo” V ° / 


cd 


Q 4 ^1850+0244 

A ^ O / 

10 '^ ‘a 

lo" 


Companion: 

o He WD 

□ COorONeMgWD 
A Ultra-light 


0.0 0.1 0.2 0.3 0.4 0.1 

Median Companion Mass (Solar masses) 


Figure 6. Left: Eccentricity as a function of orbital period for Galactic binary pulsars from the ATNF pulsar catalog. The solid line shows 
the relation from the HeWD-MSP binary formation model of IPhinnev fc Kulkarnil d 19941) and the dotted lines show the 95% confidence 
intervals for that model. Crosses with error bars show the values and uncertainies at a 95% confidence level for the binary systems in 
this work. As the measured eccentricity for PSR J1902H-0300 is consistent with zero it is shown as an upper limit. Right: Companion 
mass versus orbita l perio d for Galactic field binary pulsars from the ATNF pulsar catalog. The dotted line shows the relation predicted by 
ITauris &; SavonTi^ d 19991) for MSP-WD systems formed from LMXBs. Crosses show the four binary systems presented in this work. Other 
symbols represent the type of companion as tabulated in the ATNF pulsar catalog. 


MSPs and the addition of the PALFA discovered MSPs 
represents a 70% increase in the number of MSPs with 
DM> lOOpccm”^. Fourteen of the PALFA-discovered 
MSPs have DM/P ratios larger than 30pccm“^ ms“^, 
compared to 16 non-PALFA Galactic MSPs. PALFA dis¬ 
coveries have thus almost doubled the number of known 
Galactic MSPs with DM/P ratios larger than this value 
(Figured]). Here we discuss some of the notable proper¬ 
ties of each of the five MSPs. 

3.1. PSR J0557+1551 

PSR J0557-I-1551 is the most rapidly rotating of the 
five MSPs with a period of 2.55 ms. It is in a 4.9-day 
binary orbit and it has a DM of 103pccm“^. Its DM 
implies a distance, estimated f rom the NE2001 mode l of 
Galactic free-electron content (|Gordes fc Laziol[2002ll . of 
2.9 kpc. The minimum companion mass (i.e. the mass 
implied by the measured mass function assuming an or¬ 
bital inclination angle of z = 90° and a pulsar mass of 
1.4 Mq) is 0.20 Mq and the median companion mass (i.e. 
i = 60°) is 0.23 Mq. Our timing solution has an rms 
residual of 4.9/rs and spans 2.4years. PSR J0557-I-1551 
is located in the Galactic anticenter and so is in a region 
with a relatively low density of known pulsars. Since its 
rms residuals and TO A errors a re close to those in ex - 
isting pulsar timing arrays (e.g. IDemorest et’aTI I2013I1 . 
along with the fact that it is located in an area where 
pulsar timing arrays pres ently have low sensitivity (e.g. 
[Arzoumanian etaLl 120141) . it may be a good candidate 
for use in gravitational wave detection. 

3.2. PSR J1850+02U 

PSR J1850-1-0244 is the highest DM MSP yet discov¬ 
ered (540pccm“^); its DM-derived distance is 10.4kpc, 
the furthest for any known Galactic MSP. It has a spin 
period of 4.48 ms and is in a 0.74-day orbit with a com¬ 
panion with minimum and median masses of 0.07 Mq 
and 0.09 Mq, respectively. The timing solution for this 
pulsar has an rms residual of 57 /rs and spans 10 months. 


The profile of PSR J1850-I-0244 shows a clear scatter¬ 
ing tail. It also shows a slow rise implying that the under¬ 
lying pre-scattered profile is quite broad. We fit the pro¬ 
file with a wrapped exponentially modified Gaussian, i.e. 
an underlying Gaussian profile convolved with an expo¬ 
nential decay, caused by the scattering, that is wrapped 
at the period of the pulsar. We show the pulse profile in 
Figure [5| along with the resulting best-fit exponentially 
modified Gaussian model. The best-fit exponential de¬ 
cay timescale was found to be 1.02±0.04ms; this is close 
to the value predicted for scattering at 1.4 GHz by the 
NE2001 model of 1.3 ms. The best-fit width parame¬ 
ter of the underlying Gaussian is cr = 0.41 ± 0.02 ms 
which corresponds to a width at 50% of the peak of 
W50 = 0.97ms (a ^ 20% duty cycle). We caution 
that the exponentially-convolved Gaussian fitting func¬ 
tion is only an ap proximation to th e scattering of the ac¬ 
tual pulse shape (|Bhat et Id] 120041) . and that the above 
quoted uncertainties in the scattering time and intrinsic 
pulse width are solely statistical, and do not take into ac¬ 
count the likely dominant systematic uncertainties from 
such an approximation. 


3.3. PSR JI 9 O 5 +O 453 


PSR J1905-I-0453 is an isolated pulsar with a spin pe¬ 
riod of 6.09ms, a DM of 183pccm“^, and a DM-implied 
distance of 4.7 kpc. Its timing solution has an rms resid¬ 
ual of 33 IIS and spans 2.1 years. 

PSR J1905-I-0453 has the second longest period out of 
all known isolated MSPs (PSR J1730—2304, with a spin 
period of 8.12 ms has the longest period among isolated 


MSPs: [Lorimer et al.lll99^ . It also has the lowest E out 
of any known isolated MSP (1 x 10^^ergs“^). Isolated 
recycled pulsars are thoug ht to originate from d isruption 
due to a second super nova (iLorimer et al.ll2004f) . ejection 
from a triple system (|Freire et al.ll201~), or ablation of 
the companion by the pulsar wind ( Fruchter et al.l[l990t 
iBailes et al.ll20FUl . The short spin period requires a com- 










































panion lifetime that is much longer than that of stars 
massive enough to explode as supernovae. Ablation re¬ 
quires enough spin-down luminosity to power a pulsar 
wind that can disintegrate its companion. The question 
is then whether or not the low E of PSR J1905-I-0453 
would have been sufficiently large to power ablation ear¬ 
lier in its life as an MSP. 

The P, and thus E, of PSR J1905-I-0453 will be bi¬ 
ased by the accele ration in the Galactic plane. Using 
the formulae from iNice fc Tavlod (|1995f) . we find that 
the motion of PSR J1905-I-0453 in the Galactic plane 
causes an effective P of 1.3 x 10“^^ and that the ac¬ 
celeration in the direction perpendicular to the plane is 
negligible. The corrected P and E are thus 4.0 x 10“^^ 
and 7.0 x lO^^ergs”^. 

The P of PSR J1905-I-0453 could also be influenced 
by the Shklovskii effect whereby the space velocity 
relative to the observer causes an apparent contribu¬ 
tion to P. Since the proper motion measurement for 
PSR J1905-I-0453 is only an upper limit, the contribution 
from the Shklovskii effect cannot be calculated. However, 
if we assume that the observed P is caused exclusively 
by the Shklovskii effect, it implies a transverse velocity 
of 200km/s. This is a reasonable velocity; according to 
the ATNF pulsar catalog, ^ 40% of pulsars with mea¬ 
sured proper motions have inferred transverse velocities 
equal or greater than 200km/s. The true P, and hence 
P, of PSR J1905-I-0453 could therefore be much lower 
than the value in Table [T] 

For now, let us assume that the contribution from the 
Shklovskii effect is minimal. The spin-down luminosity 
of PSR J1905-I-0453 was of course higher in the past. 
If we assume pure magnetic-dipole braking, a constant 
magnetic held, and that the MSP was recycled at most 
14Gyr ago (the age of the Universe), the fastest period 
and highest E that PSR J1905-I-0453 could have had 
at formation would be 3.34ms and 8 x 10^^ergs“^, re¬ 
spectively. This spin-down luminosity is comparable to 
those of several black widow pulsars, the least energetic 
of which have P ^ 1 x lO^'^ergs"^ (jRobertsll2013[) . Since 
black widow pulsars are thought to be in the process of 
ablating their companion, this suggests that it may have 
been possible for PSR J1905+0453 to have done so given 
its spin-down power if it originated in a binary similar 
to those of known black widow pulsars and the contri¬ 
bution of the Shklovskii effect to its measured P is not 
dominant. 

3.4. PSR J1902+0300 

PSR J1902-I-0300 has a period of 7.80 ms, a DM of 
254pccm“^, and a NE2001 implied distance of 5.9kpc. 
It is in a 2.4-day orbit and its companion has a minimum 
mass of O.I 4 M 0 and a median mass of O.I 6 M 0 . Our 
timing solution has an rms residual of 14 /xs and spans 
2.3 years. 

3.5. PSR JI 943 + 22 IO 

PSR J1943-I-2210 has a period of 8.31ms, a DM of 
174pccm“^, and a NE2001 implied distance of 6.2kpc. 
It is in a 8.3-day orbit with a minimum companion mass 
of 0.28 Mq and a medium mass of 0.34 Mq. Its timing 
solution has an rms residual of 13 /iS and spans 1.2 years. 


3.6. Counterparts in other wavelengths 

Because the pulsars presented here are quite faint and 
distant we would probably not expect counterparts in 
other wavelengths to be easily detected. The brightest 
white dwarfs have opti cal absolute magnitudes of ~ 10 
(|Garroll fc Ostliell2007f) . The closest of the five binary 
systems in this work, according to their DM-derived dis¬ 
tances, is PSR J0557+1551 at 3 kpc. At this distance a 
Mv = 10 white dwarf would have an apparent magni¬ 
tude of ^ 22. This is quite faint and it is likely that the 
actual companion is fainter than this. 

To see whether we expect to see gamma-ray counter¬ 
parts for the MSPs, we can compare the spin-down en¬ 
ergy incident at the Earth, P/of the pulsars here with 
those of the Fermi detected MSPs. According to the 
ATNF pulsar catalog, most Fermi-detected MSPs have 
P/P2 > 10^3 ergs-^kpc-2. Only PSR J0557-bl551 has 
a E/D"^ higher than this, and it is possible that a signif- 
cant portion of the measured P is due to the Shklovskii 
effect. Furthermore, the other four MSPs are located 
in the inner Galactic plane where the gamma-ray back¬ 
ground is high. 

Indeed, we find no coincident optical, infrared, X- 
ray or 7 -ray sources in the SIMBAD^'^ or HEASARC^® 
datab ases or in the Ferm i LAT 2-year Point Source Gat- 
alog (jNolan et al.ll20l'^ for any of the MSPs reported 
here. 


3.7. Nature of the Binary Companions 

None of the four binaries presented here show any evi¬ 
dence for eclipses or matter in the system in their timing 
residuals, which suggests that none of the companions 
are likely close to filling their Roche lobes. Given the 
orbital separation and masses of the components of a bi¬ 
nary system, the s ize of the Roche lobe of the com panion 
can be calculated (jTauris fc van den Heuve]|[2006D . If we 
assume a typical inclination of 60° and a pulsar mass of 
1.4 M 0 , the sizes of the companion Roche lobes for all 
four binary systems can accomodate both white dwarfs 
and main sequence stars. 

In Figure [ 6 ] we plot the binary parameters of the four 
binary MSPs presented here along with those from all 
Galactic-field binary MSPs, taken from the ATNF pul¬ 
sar catalog. The plots of orbital period versus eccentric¬ 
ity and companion mass versus orbital period allow us to 
com pare the binary param e ters t o the predicted relations 
from l Phinnev fc Kulkarnil (jl994f l and iTauris fc Savoniid 
([19^ for helium white-dwarf (HeWD)-MSP binaries 
that evolved through long-term, stable, mass transfer 
from a red giant progenitor onto a neutron star in a low- 
mass X-ray binary (LMXB). 

The median companion masses implied by the mass 
functions of PSR J1902-b0300 and PSR J0557+1551 
are in the range expected for HeWD companions 
(jTauris et al.ll2ni^ . They also fall on the orbital period- 
eccent ricity relation predicted by iPhinnev fc Kulkarnil 
(Il994f) and near th e companion mass - orbita l period re¬ 
lation predicted bv lTauris fc Savoniid (1199911 . Thus, the 
likely nature of their companions are HeWDs. 

http: //simbad.u-strasbg. fr/slmbad/ 
http://heasarc.gsfc.nasa.gov/ 
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The other two binary systems fall slightly outside of the 
parameters expected for a system with a HeWD compan¬ 
ion and so formation scenarios other than LMXB evolu¬ 
tion into a HeWD may have occured. PSR J1943-I-2210 
has a companion mass that is higher than expected 
for a HeWD and f alls o utside of the relation from 
iPhinnev fc Kulkarnil (|1994f) with a lower than expected 
eccentricity. The high companion mass implies that its 
companion could be a carbon-oxygen white dwarf. It 
would therefore be c onsidered an intermediate-mass bi¬ 
nary pulsar lIM BP: iC amilol I 1996H a class of currently 
about 20 objects (|Tauris et al.ll201^ . 

PSR J1850-1-0244, on the other hand, has a lower mass 
than expected for a HeWD and a higher than predicted 
eccentricity. In the Me — Pb plot of Figure El it falls 
close to the systems classified as “Ult ra-light”, man y 
of which are considered black widows (IRobertsi 1201311 . 
raising the possibility that it has interacted with its 
companion through ablation and accretion processes, as 
the black-widow systems are thought to have done. Of 
course, the true companion mass of PSR J1850-I-0244 
could be higher if the system has a low inclination. 
The relativistic MSP-HeWD binaries PSRs J0348-I-0432 
and J1738-I-0333, which like PSR J1850-f0244 have or¬ 
bital periods less than a day and median compan¬ 
ion masses ^ 0.1 Mq, have measured true companion 
masses determined from their post-Keplerian parame¬ 
ters and white dwarf spectroscopy. PSR J0348-I-0432 
has a companion mass of 0.1 72 Mq and an inclination 
of 40° (|Antoniadis et al.l 1201311 . PSR J1738-I-0333 has a 
companion mass of 0.18 1 Mq and an inclination of 32° 
(|Antoniadis et all 1201211 . These com panion masses are 
close t o those predicted for HeWDs bv lTauris fc Savoniiel 
(|1999H . If we assume that PSR J1850-I-0244 has a HeWD 
companion with a mas s of 0.19 Mq as predicted by 
iTauris fc Savonii^ (11999( 1. the inclination of the binary 
system would be 24° for a pulsar mass of 1.4 Mq and 
30° for a pulsar mass of 2.0 Mq. This is not an unrea¬ 
sonably fortuitous orientation and so it is plausible that 
PSR J1850-1-0244 has a low inclination and a standard 
HeWD companion. 

3.8. Future prospects for mass measurements 

Masses of the components in a MSP bin ary can be mea¬ 
sured using the Shapiro delay effect ('e.g.lRyba fc Tavloii 
119911: iSolaver et al.l (200^ iDemorest et al.1 12010( 1. The 
Shapiro delay is parameterized in timing models by the 
post-Keplerian parameters r and s. The Shapiro range, 
r, and shape, s, parameters are equal to the companion 
mass in units of light seconds and the sine of the inclina¬ 
tion, respectively. 

In the ELLl model, when binary systems have mea¬ 
sured ei = esinw greater than zero and £2 = ecosw 
consistent with zero, it may be the case that the eccen¬ 
tricity, e, is actaully zer o and the non-zer o ei is due to 
a Shapiro delay signal (iLange et al.ll200ill . For two of 
the four binary MSPs in this work, PSRs J0557+1551 
and J1850+0244, ei > 0 and £2 = 0. For these sys¬ 
tems, given an assumed inclination, we can calculate the 
companion mass (i.e. Shapiro r) implied by the mea¬ 
sured £1 assuming that the true ec centricity is consis - 
tent with zero (see Appendix A2 of (Lange et al.l [2001(1 . 
If the implied companion masses are in a reasonable 
range given the measured mass functions, this would 


give evidence that the measured £1 is due to Shapiro 
delay. For PSR J0557+1551, the lowest implied com¬ 
panion mass (i.e. at high inclination) is ~ 2 Mq and for 
PSR J1850+0244 it is ^ 2.5 Mq. Since the measured 
mass functions for these two MSPs would require unrea¬ 
sonably high pulsar masses in order for their companions 
to be so heavy, this implies that the measured values of 
£1 for PSRs J0557-I-1551 and J1850-I-0244 are not domi¬ 
nated by Shapiro delay signals. 

In order to test for a Shapiro delay signal in each of the 
binary puls ars, a timing model was fit using the TEMPO 
DD model (|Damour fc Deruelle( I1985L (1986(1 which in¬ 
cludes post-Keplerian parameters in the timing model. 
We fit the TOAs using TEMPO with the best-fit param¬ 
eters listed in Table [1] with the addition of a test Me 
and cos*. We searched a grid of Me and cos* ranging 
from 0.0 < Me < 1.0 and 0.0 < cos* < 1.0. For all but 
PSR J0557-I-1551, all points in the grid provided accept¬ 
able fits at greater than a 95.4% level and so no limits 
could be placed. For PSR J0557-I-1551, the left panel of 
Figure [7] shows the upper-limit contours placed on M 2 
and cos*. It is clear that a small portion of reasonable 
binary system parameters are inconsistent with the data 
at a Icr confidence level, namely at high inclination and 
high pulsar mass. 

These current limits are not very constraining, but fu¬ 
ture observations may yield interesting results. To inves¬ 
tigate this, we constructed a simulated set of TOAs with 
uncertainties and rms residuals of approximately 1 /iS. 
The simulated timing campaign spanned a year with a 
one month cadence and included ten 2.5-h observations 
that each occured at one of the ten conjunctions of the 
system that will be visible at Arecibo in that time span. 
The right panel of Figure 0 shows the limits that could 
be placed with such a campaign given a pulsar mass of 
1.4Mq and an inclination of 80° which imply a com¬ 
panion mass of 0.20Mq. Note that if PSR J0557+1551 
has a more favorable binary configuration (i.e. higher 
component masses or inclination) or if the timing pre¬ 
cision could be improved to better than l-/is then more 
stringent measurements could be made with a campaign 
similar to that simulated. 

4. CONCLUSIONS 

The search for MSPs is heavily motivated by the 
unique and exciting studies that are enabled by finding 
more systems that provide laboratories for such science. 
PALFA provides a unique window on the MSP popula¬ 
tion, discovering more distant and fainter MSPs than any 
other past or current survey. We have presented here the 
discovery and follow-up timing of five MSPs discovered 
in the PALFA survey as well as discussed their interest¬ 
ing properties. One MSP, PSR J1905+0453 is isolated 
and has the lowest spin-down luminosity yet measured 
for an isolated MSP. Three of the binary MSPs likely 
have HeWD companions and one, PSR J1943-T2210, may 
have a carbon-oxygen WD companion. At least one of 
the five, PSR J0557-T1551, has the potential to allow 
mass measurements and may be useful in pulsar timing 
arrays. PSR J1850+0244 is the highest DM MSP yet 
discovered with a DM of 540. Further discoveries in the 
ongoing PALFA survey will expand our knowledge of the 
MSP population deep into the Galactic plane as well as 
help constrain fundamental physics and astrophysics. 
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Figure 7. Maps of companion-mass and inclination angle for PSR J0557-I-1551. The left plot shows binary configurations allowed by 
current timing data and the right shows configurations for a simulated data set with Mp = 1.35Mq and i = 80° (see text). The black cross 
denotes the simulated binary configuration. The black contours enclose 68.3% (Ifr), 95.4% {2a), and 99.7% (3cr) of allowed configurations. 
The red dashed lines denote the pulsar mass for the configuration. 
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Table 1 

Timing Parameters for five MSPs 


PSR J 

0557-1-1551 

1850-1-0244 

1902+0300 

1905+0453 

1943+2210 

Measured Parameters 

Right ascension (J2000) 

05 : 57 : 31.44918(9) 

18 : 50 : 41.683(1) 

19 : 01 : 59.6168(1) 

19 : 04 : 59.3848(2) 

19 : 43 ; 16.4910(1) 

Declination (J2000) 

15 : 50 : 06.046(8) 

02 : 42 : 57.44(3) 

03 : 00 : 23.324(5) 

04 : 51 : 54.953(8) 

22 : 10 : 23.151(2) 

Total proper motion, fiT (masyr~^)“ 

< 80 

- 

< 77 

< 96 

- 

Spin frequency, u (s~^) 

391.18012787837(8) 

223.221059721(3) 

128.25813729784(5) 

164.1397912192(1) 

196.6887956594(2) 

Frequency derivative, z> (s~^) 

-1.125(3) X 10-1® 

-8.1(2) X 10-1® 

-7.52(2) X 10-1® 

-1.54(6) X 10-1° 

-3.4(1) X 10-1° 

Dispersion measure, DM (pccm“^) 

102.5723(5) 

540.068(9) 

253.887(2) 

182.705(4) 

174.089(3) 

Orbital period, (d) 

4.846550440(4) 

0.73543015(1) 

2.399218509(6) 

- 

8.31148080(2) 

Projected semimajor axis, x = asini/c (s) 

4.0544597(8) 

0.452355(9) 

1.842553(1) 

- 

8.064088(2) 

Time of ascending node, Tasc (MJD) 

56349.2459263(2) 

56468.179498(6) 

55989.555577(1) 

- 

56392.7364900(7) 

ei = esincj 

9.3(4) X 10-® 

9(3) X 10-® 

3(2) X 10-® 

- 

2.5(8) X 10-° 

£2 = 6 cos OJ 

-9(60) X 10-® 

-4(40) X 10-® 

4(10) X 10-1 

- 

-1.4(4) X 10-° 

Derived Parameters 

Galactic longitude, 1 (degrees) 

192.68 

35.26 

36.81 

38.81 

58.44 

Galactic latitude, b (degrees) 

-4.31 

1.40 

-0.97 

-0.79 

-0.73 

Period, P (ms) 

2.5563670767830(5) 

4.47986404710(6) 

7.796776259723(3) 

6.092367929629(5) 

5.084173689952(5) 

Period derivative, P 

7.35(2) X 10-21 

1.62(4) X 10-1® 

4.57(1) X 10-®° 

5.7(2) X 10-®i 

8.7(3) X 10-®i 

Eccentricity, e 

9.3(4) X 10-® 

9(3) X 10-® 

3(2) X 10-° 

- 

2.9(7) X 10-° 

Longitude of periastron, oj (deg) 

91(4) 

90(20) 

80(30) 

- 

120(10) X 10® 

Time of periastron passage. To (MJD) 

56350.47(5) 

56468.37(5) 

55990.1(2) 

- 

56395.5(2) 

Mass function (Mq) 

3.046605(2) x 10-® 

1.8375(1) X 10-^ 

1.166818(2) X 10-® 

- 

8.150639(6) X 10-® 

Minimum companion mass (Mq)^ 

0.20 

0.07 

0.14 

- 

0.28 

Median companion mass (Mq)*^ 

0.23 

0.09 

0.16 

- 

0.34 

Surface dipolar magnetic field, B (G) 

1.388(2) X 10® 

8.6(1) X 10® 

6.040(8) X 10® 

1.89(3) X 10® 

2.13(3) X 10® 

Spin-down luminosity, E (ergs~^) 

1.738(5) X 10®-! 

7.1(2) X lO®'! 

3.81(1) X 10®® 

1.00(4) X 10®® 

2.61(8) X 10®® 

Characteristic age, Tc (yr) 

5.51(2) X 10® 

4.4(1) X 10® 

2.705(7) X 10® 

1.69(6) X 10®° 

9.3(3) X 10® 

DM distance, D (kpc)*^ 

2.9 

10.4 

5.9 

4.7 

6.2 

Transverse velocity, Vt (kms~^)°' 

< 1100 

- 

< 2100 

< 2100 

- 

Observation Parameters 

1400 MHz mean fiux density, 61400 (riiJy) 

0.050(6) 

- 

0.113(4) 

0.117(9) 

0.04(1) 

1400MHz mean radio luminosity, I /1400 (iriJy kpc^)® 

0.4 

- 

~ 3.9 

~ 2.6 

~ 1.5 

Discovery observation date (MJD) 

55907 

56466 

55863 

55817 

56389 

Timing epoch (MJD) 

56346.00 

56554.29 

56242.00 

56242.00 

56404.00 

Start epoch (MJD) 

55963.08 

56546.71 

55989.20 

56067.92 

56391.07 

Finish epoch (MJD) 

56841.65 

56856.14 

56828.22 

56831.16 

56828.25 

Number of TOAs 

162 

160 

270 

278 

147 

Typical Integration time per Arecibo TOA (min) 

10 

2 

2 

2 

2 

Typical TOA uncertainty (/rs) (Jodrell/AO-search/AO-fold) 

23/7/3 

130/50/- 

80/15/15 

160/30/25 

60/10/10 

TOA error scale factor, EFAC (Jodrell/AO-search/AO-fold) 

1.58/1.95/1.67 

1.67/1.13/- 

1.19/1.06/1.07 

2.34/1.01/0.98 

2.19/0.91/0.89 

rms post-fit residuals (/r.s) 

4.87 

57.12 

13.92 

32.92 

13.31 


Note. — Numbers in parentheses are TEMPO reported Icr uncertainties. 
^ 95.4% (2(j) upper limits. 

^ Assuming a pulsar mass of 1.4 Mq and an inclination of 90°. 

° Assuming a pulsar mass of 1.4 Mq and an inclination of 60°. 

As implied by the lCordes & Laziol 11200211 NE2001 DM-distance model. 
® 1/1400 = 5'i400 















